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Abstract 



H ' We demonstrate that CP— violation in the Majorana mass matri- 

ces of the heavy neutrinos can generate a CP— asymmetric universe. 
The subsequent decay of the Majorana particles generates a lepton 
number asymmetry. During the electroweak phase transition the lep- 
ton asymmetry is converted into a baryon asymmetry which survives 
down to this time. 
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The cosmological baryon excess can be generated from the initial con- 
dition B = through the baryon number violating decays of very heavy 
particles [|TJ. This requires violation of both C and CP and, in addition, 
these processes must happen during an epoch which brings the universe out 
of equilibrium. This mechanism has two major drawbacks. In most grand 
unified theories, baryon changing decays conserve (B — L) and anomalous 
electroweak interaction induced by sphalerons can wash out the asymmetry 
[§]. Second, the asymmetry generated in most models is much smaller than 
the ratio 
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required by observations and big-bang cosmology. A second approach relies 
on the violation of baryon number through the global (B — L)— anomaly 
which together with C— and CP— violation generates a baryon excess at the 
electroweak phase transition, provided it is first order ||. The search for a 
consistent extension of the standard model is under active investigation. 

At present, a very attractive scenario consists of generating a lepton 
number asymmetry through lepton number violating decays of heavy neutri- 
nos. The baryon asymmetry is then generated from the lepton asymmetry 
during the electroweak phase transition, mentioned above. Some time ago, 
Fukugita and Yanagida proposed this mechanism [4-7] where the lepton num- 
ber violation was introduced explicitly through a Majorana mass term of the 
right-handed neutrinos. In all these models the CP— violation is introduced 
through the interference of tree-level and one-loop diagrams in the decays of 
the heavy neutrinos. 

Attempts were made |8|, || |10| to generate the baryon asymmetry in mod- 
els, where the CP— violation comes through the mixing of scalars. The pos- 



sibility of mixing scalar particles in order to obtain large CP— violation was 
discussed some time back ||. A similar idea was recently exploited in the 
context of GUTs |J |H]], where the oscillation of scalar particles was proposed. 
In these scenarios, the GUTs were extended to include new scalars which al- 
low such interactions and introduce e-type CP— violating effects in the mass 
matrix of the scalars. In the above models CP— and baryon-violation oc- 



cur in the same diagrams. In a recent article ||11| , e-type effects have been 
considered in a different way. Through the scattering of heavy bosons it has 
been shown that it is possible to create first a CP— asymmetric universe (to 
be precise CP— asymmetric density of scattering states) before the heavy 
particles decay. When they begin to decay they generate the baryon asym- 
metry. The main difference of this mechanism from the previous ones lies in 



the origin of CP— violation. In the latter mechanism [O] the CP— violation 
arises from the phases in the couplings of the matter-anti-matter oscillations 
in which a quantum number changes by two units, in contrast to references 



[^, [| [TD[ where the CP— violation occurs in decays of the heavy particles and 
the quantum numbers change by one unit. If this oscillation violates baryon 
number by two units, and at the same time there is CP— violation, then 
the oscillation will generate a CP— asymmetric and also baryon asymmetric 
universe. 

In the case of heavy Majorana neutrinos, there exists one loop self-energy 
diagrams introducing a correction to the Majorana mass of the particles. 
These diagrams will not conserve CP and will produce neutrino-anti-neutri- 
no oscillations violating lepton number by two units, which can contribute 
to the lepton number asymmetry of the universe. In the case of three gen- 
eration models, the CP— violating phases which enter in the self energy di- 
agrams are usually not the ones appearing in the interference terms of the 



decays of the heavy neutrinos | I2| . Thus this contribution of C— , CP— and 
(B — L)— violation to leptogenesis is distinct from the ones considered in the 
literature. In addition, the contributions are larger than the earlier ones for 
several choices of parameters. 

We consider a theory with the same particle content as in ref @], that is, 
the particle content of the standard model plus three Majorana neutrinos Ni 
with i = l,2, 3, one for each generation. The new Majorana neutrinos are 
singlets under the standard model, i.e. SU(3) C x SU{2) L x U{l)y- The new 
interactions of these neutrinos with leptons and the Majorana mass term are 
given by, 

C = MijNfNj + haNilaft + h.c. (1) 

where <fi is the standard model higgs doublet, which breaks the electroweak 
symmetry and gives mass to fermions. 

The lepton number and (B — L) are violated through the Majorana mass 
of the right-handed neutrinos iVVs. In all other interactions lepton number 
is conserved. We shall work in the density matrix formalism |13| to estimate 



the contribution of neutrino-anti-neutrino oscillations to the lepton number 
generation of the universe. The neutrino state \N > and the anti- neutrino 
state \N C > can be distinguished by their decay properties. The state \N > 
decays into a light lepton U and a higgs <j)\ whereas the state \N C > decays 
into a light anti-lepton l\ and a higgs 0. The only interaction which re- 
lates the two states and violates lepton number is the Majorana mass term, 
which also causes mixing of the two states. We shall discuss the question of 
CP— violation in such a mixing later on. We shall assume CPT— invariance, 
but violation of CP. For the sake of simplicity, we consider two generations 
of Majorana neutrinos where the indices i and j take the values 1 and 2. We 



define a basis of states as (iVf N 2 N\ N 2 ) and the Hamiltonian takes the 

form, 

( H n # 12 \ 
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Without loss of generality, we work in the basis where the Majorana mass 
matrix is diagonal and real. Then the terms Hu and H 22 are the mass 
terms in equation (1). In this basis, the element H 12 is generated through 
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the diagrams of figure 1. Similar diagrams produce H 2i = M ; 

m 21 = m? 2 , r 2 i = rj 2 . 

We define the 4x4 density matrix in the basis of states defined above. 
It satisfies the equation, 

dp 



dt 
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The solutions are obtained by diagonalizing H. The expectation value of any 
observable, 

\N{ >< iV 2 c | . 



Of =< f\ 
is given by, 
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where |/ > is the final states of leptons or anti-leptons produced through the 
decays of these density of state \N >— > l<jy and |iV c >^ l c <p. The observables 
corresponding to leptons and antileptons in the final states are, respectively, 

\ 
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and Or oc 



fh al h* al h al h* a2 0\ 
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Hence a lepton asymmetry is generated through the CP- violating phase and 
is given by, 

S = ^— ?• (3) 

The above asymmetry S calculated later on is different than the asymme- 



try generated from the decay of heavy neutrinos |]T0[ , because the elements 
of the density matrix already contain CP-violation. In the terminology of 
K— meson decays it originates in the mass-matrix, that is the density matrix 
described above (similar to that of ref.[IT|)Q. 



In the present case, the physical states are pure states and can also be 
described in the wave function formalism. For the mass matrix under con- 
sideration, we give the mass eigenstate for one of the eigenvalues as, 

V>i « |JVi + N? > +b\N 2 > +d\N' > 

The decay widths for the decays of ip\ into leptons and antileptons is given 
by 

Ti oc \h al + b h a2 \ 2 
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2 We thank Dr. H. So for pointing out that the specific mechanism in rcf. lllj is sup- 
pressed by the general property of GUTs, according to which (B — L)-changing operators 
are suppressed relative to the (B + L)-changing operator. In ref. [ plj the suppression is 
implied by the conservation of weak isospin. 



derived from the Hamiltonian perturbatively. In a similar way we can now 
write down the decay rates for the states ip%, ^3 and 1JJ4. The quantities b and 
d are determined form the loop diagram in figure 1. The final asymmetry 
obtained either in the density matrix or the wave function formalism. The 
oscillation of a j = 2 neutrino to an % = 1 neutrino and its decay to the ath 
light lepton or antilepton produces the asymmetry, 

We shall include this lepton asymmetry in the evolution of the lepton number 
determined by the Boltzmann equation. 

The number density of the left-handed leptons n\ will evolve in time 
following the equation, 



dni I iu 

— + 3H m 

dt 



A l n j T i + -(l + e')n l T t -V i 



(5) 



The second term on the left side of the equation comes from the expansion 
of the universe, where H is Hubble's constant. We now discuss the origin of 
the various terms on the right of this expression. 

The first term contains the number density rij of initial states \Nj >, 
which is converted to the final state |iV? >, which subsequently decays into 

h 2 

light leptons with a decay width I\ = j^Mi. Here Mi is the Majorana mass 
eigenstate of the i-th neutrino. The factor A; contains the CP— conserving 
and CP— violating terms of the \N > —\N C > oscillations. The second term 
in the expression ([5]) arises from the CP— violation in the decays of the heavy 
leptons. The interference of the tree level diagram and the one loop penguin- 
type diagrams in the decays of the neutrinos iV, contains the CP— violating 
phase and the absorptive part of the integral. This contribution is the usual 
one which was studied in the literature in detail [f|, |5|, ||. The measure of the 
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CP— violation is now given by e'. In the absence of the new contribution, i.e., 
of the term proportional to Aj, the asymptotic value of the lepton number 
asymmetry is proportional to e'/g*. 

The evolution of the light anti-leptons will be given by similar equations, 
with A/ — > Aj, e 1 — > —e' and / — » I. This way we arrive at a Boltzmann 
equation for the lepton difference n^s = ni — nj, where s = #*7i 7 is the entropy 
density, g* the number of active degrees of freedom and n 7 is the number 
of photons. In the difference occurs a term proportional to 5 = A/ — Aj 
computed in equation (^). 

There are also various processes which deplete the number densities of the 
left-handed leptons. They come from the reactions I + <p + — > l c + 0. We can 
denote by Di the depletion of lepton and Dj the depletion of the antileptons. 
The difference which occurs in the asymmetry equation is, 

V = D l -D T = (25 + eV^r* + a v ) + ^{T, + a~ v )— (6) 

n 7 

The term with Tj comes from the \N > or |iV c > neutrinos in the intermedi- 
ate state and d v is produced from the remaining intermediate states. Since 
this contribution is, in general, smaller than F it will be neglected. Direct 
and indirect CP— violation in the intermediate states produce the first term. 

The difference \A(l(j) c -\ > l c 4>)\ 2 - \A(l c (f)-\ ► l(j) c )\ 2 will also contribute to 

the depletion through the term proportional to M '~ M ' = n '~ n ' ■ When d v is 
neglected, the last term in this expression is due to the recombination of the 
left-handed leptons and the higgs to form the heavy neutrinos. All these 
processes are rapid enough to occur in thermal equilibrium and hence the 
number density entering here is the equilibrium number density n eq of N. 

We now introduce Yj = rii/s, Y eq = n eq /s, the dimensionless variable, 
x = Mi/T and make use of the relations, t = 1/2H = x 2 /2H(x = 1) with 



M 2 

H(x = 1) = 1.73y/g^-jrr~ (from now on we assume that Mi < M2 < M3 and 



-M Pl 

the lepton number asymmetry is generated just before T ~ Mi). We also 
define the parameter, 

_ r,(x = 1) 

which is the measure of the out-of-equilibrium condition. For the usual ex- 
tensions of the standard model, one can take g* = 400. With these redefini- 
tions we obtain a simplified equation for the evolution of the lepton number 
asymmetry of the universe, 

^ = Kx 2 [{25 + e') (Y t -Y^)- Y^g*n L ] . (7) 

The Yi satisfy the Boltzmann equations, 

c\Y 

^l = -Kx 2 (Y l -Y e i). (8) 

To solve these equation we assume that although the universe is out-of- 
equilibrium, it stays close to it, so that ' *T — - = 0. Then it is easier to 
solve the equations and obtain an estimate of the lepton number generated. 
For very large time, the solution of a differential equation has an asymptotic 
value, which is almost independent of the constant K for 1 < K < 50, and 
is given approximately by, 

til = -(26+ J). (9) 

9* 

If 5 = 0, we recover the expression for the scenarios 0, |5|], where the heavy 
neutrinos can only decay. We compare next the two CP— violating contribu- 
tions in these scenarios. 

The CP— violation in the decays of heavy neutrinos has been discussed 
in the literature and for the model under consideration it is given by [|], |5], §], 

= 1 \m(h al h aj h* 01 h* pj ) Mi 

2vr h al h* al Mj' [ ' 



From this expression it is clear that the present contribution 5 is different 
from the contribution considered earlier, given be e'. There exists various 
choices of parameters for which 5 ^> e'. In all these cases, the matter-anti- 
matter oscillations will first create a CP— asymmetric universe with unequal 
densities of Majorana neutrinos and antineutrinos, whose decays generate a 
lepton asymmetric universe. 

In the presence of sphalerons, this lepton number asymmetry will be 
converted to the baryon number asymmetry. After the electroweak phase 
transition, the baryon number of the universe will be given by, 

/ 8N g + AN H \ 1 25 

HB = ~ [22N g + l3N H ) "<*"'> ~ 3 nL ~ %? 

where Nh is the number of Higgs doublets and N g the number of generations. 

For h\z ~ hu and g* ~ 400, the required amount of baryon asymmetry can 

be generated with a CP— violating phase 5 ~ 10~ 5 — 10~ 7 , which is a very 

natural choice. 

We have shown that CP— violation in the mass matrix of heavy Majorana 

neutrinos produces a CP and lepton asymmetric universe. The same mass 

terms break the (B — L) symmetry and produce a net excess of this quantum 

number which survives all the way down to our epoch. The asymmetries 

created on the mass matrix are relatively large and are transformed into an 

excess of Baryons, when the electroweak phase transition takes place. 
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Figure 1: One loop diagrams contributing to the CP violating |JVi >- 
\N% > transition. 
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